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Abstract 
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A  numerical  method  is  presented  for  determining  the 
bending  stresses  in  elastic  cantilever  plates.   Southwell's 
formulation  of  the  classical  small  deflection  plate  theory 
is  employed.   The  method  is  applied  to  obtain  numerical 
solutions  for  rectangular  and  swept  back  plates  subjected 
to  a  variety  of  surface  and  edge  loads.   Boundary  layer 
effects  are  observed  for  the  uniformly  loaded  rectangular 
plates  as  the  length  of  the  plate  decreases.   The  effects 
of  a  load  concentration  and  sweep  back  are  observed  and 
summarized  in  the  form  of  stress  concentration  factors. 
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1.  Introduction. 

The  determination  of  the  stresses  in  elastic  cantilever 

plates  is  required  for  a  variety  of  structural  applications. 

The  plate  is  clamped  on  a  segment  C  of  the  boundary  curve  B.   On 

B*  =  B  -  C,  the  normal  bending  moment  M  and  the  shear  reaction 
—     '  a  n 

V  are  specified.   The  face  of  the  plate  is  subjected  to  a  normal 

pressure  p(x,y). 

The  applied  forces  M  ,  V  and  p  and  the  thickness  of  the 

n   n 

plate  are  sufficiently  small  so  that  the  classical  small  deflec- 
tion, Lagrange -Kirchhof f  theory  of  plates  [1]  is  valid.   In  its 
conventional  formulation,  this  theory  requires  the  solution  of 
a  biharmonic  boundary  value  problem  for  the  normal  displacement 
w(x,y)  of  the  midplane  of  the  plate.   The  boundary  conditions  are 
that  w  and  its  first  normal  derivative  vanish  on  C,  and  two 
linear  combinations  of  the  second  and  third  derivatives  of  w  are 
specified  on  B' . 

Southwell's  [2]  formulation  of  the  classical  theory  yields 
a  boundary  value  problem  for  a  system  of  three  second  order 
partial  differential  equations  for  three  dependent  variables 
(j)(x,y),  u(x,y)  and  v(x,y).   The  function  <j)  is  a  solution  of  a 
boundary  value  problem  for  Poisson's  equation  and  is  determined 
independently  of  u  and  v.   The  remaining  two  equations  for  u  and 
v  are  coupled.   They  are  similar  to  the  displacement  equations 
of  equilibrium  for  either  plane  stress  or  plane  strain  [4].   The 
function  <f>  appears  in  the  inhomogeneous  terms  in  these  equations. 
The  boundary  conditions  on  B'  are  that  u  and  v  are  specified,  i.e. 
Dirichlet  conditions  are  specified  on  B'.   Thus  in  Southwell's 


he  compl  ier  derivative  boundary 

cor.     na  on  B1  of  the  classical  formulation  are  replaced  by 
nple  con :      .       is  especially  important  for  the  analysis 
paper  since  finite  difference  methods  are  easily  and 
accurately  applied  to  Dirichlet  conditions.   A  discussion  of 

"erence  and  other  approximate  methods  for  cantilever  plates 
and  references  to  previous  work  are  given  In  [1]  and  [J], 
Southwell's  formulation  has  been  extended  [5]  to  plates  of 
variable  thickness.   Numerical  solutions  are  obtained  in  [5]  by 
employing  relaxation  methods. 

In  this  paper  an  accelerated  iteration  method  is  presented 
for  solving  cantilever  problems  using  Southwell's  formulation. 
The  method  was  previously  used  in  studies  of  the  nonlinear 
deformations  of  elastic  plates  [6,7,8]  and  other  problems.   Each 
iterate,  is  determined  as  the  solution  of  a  boundary  value 
problem  for  Poisson's  equation.   The  iterates  are  numerically 
evaluated  by  approximating  this  boundary  value  problem  by 
difference  equations.   The  resulting  system  of  linear  algebraic 
equations  which  is  of  quasi-tridiagonal  form  is  solved  by  a 
direct  factoring  method  [9].   The  method  can  also  be  considered 
as  an  iterative  procedure  for  solving  the  difference  equations 
which  are  obtained  from  the  original  plate  boundary  value  problem. 
In  previous  applications  of  difference  methods  ,  coarse  meshes,  e.g. 
6  =  L/2  were  employed.   Here  6  is  the  mesh  width  and  L  Is  a 

See  e.g.  [1,3] ■ 
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minimum  diameter  of  the  plate.   The  present  computations  indicate 
that  mesh  widths  as  small  as  6  =  L/52  may  be  required  to 
accurately  estimate  the  stresses.   The  finer  meshes  are  more 
important  for  irregular  domains  such  as  the  swept  back  plate 
considered  in  this  paper;  see  Figure  1. 

In  Section  5  the  iterative  method  is  applied  to  obtain 
numerical  solutions  for  rectangular  plates  and  the  swept  back 
plate  shown  in  Figure  1  subjected  to  a  variety  of  loads.   It  is 
shown  that  the  solution  of  the  infinite  strip  cantilever  provides 
a  reasonably  good  estimate  of  the  numerically  determined  maximum 
moment  for  the  symmetrically  loaded  rectangular  plates  that  are 
analyzed  in  Section  5.   It  is  an  underestimate  and  improves  as 
the  plate  is  made  narrower.   The  infinite  strip  cantilever  gives 
an  inaccurate  approximation  for  unsymmetric  loading.   For 
equivalent  loading  conditions,  the  maximum  moment  on  the  clamped 
edge  is  larger  for  the  swept  back  plate  than  for  the  rectangular 
plate  with  the  same  aspect  ratio.   Various  stress  concentration 
factors  that  indicate  the  relative  effects  of  sweep  back  and 
unsymmetrical  loading  are  presented. 


2.  Formulation. 

The  midplane  of  the  plate  covers  the  region  R  in  the  right 
half  of  the  x,y  plane.   The  boundary  B  of  R  is  a  positively 
oriented,  simple,  closed  contour  with  arc  length  s.   The  unit 
outward  normal  to  B  is  n  =  (a,b)  and  the  unit  tangent  vector  is 
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t  =  (-b,a).   Here  a  and  b  are  the  direction  cosines  of  n.   The 

boundary  contains  the  straight  line  C:   x  =  0,  yQ  <  y  <  y^. 

The  plate  is  clamped  on  C.   On  B1  =  B  -  C,  the  normal  bending 

moment  M  =  M(s)  and  the  transverse  shear  reaction  V  =  V(s) 
n  n 

where  M  and  V  are  prescribed  functions.   An  arc  of  B'  is  free  if 
M  =  V  =.  Ci  on  that  arc.   A  pressure  p(x,y)  is  applied  normal  to 
the  x,y  plane.   The  applied  loads  p,  M,  V  and  the  region  R  are 
referred  to  as  the  data. 

The  x  and  y  bending  moments  are  denoted  by  M  and  M  and 
the  twisting  moment  by  M   =  "Myx-   They  are  related  to  the  normal 
displacement  w(x,y)  of  the  mldplane  by 


Mx  =-D(wxx+wyyh  My  =  "D(wyy  +  V  Wxx}'  Mxy  =  D(1"v)wxy'    (la'b'c) 


where  v  is  Poisson's  ratio  and  D  is  the  flexural  rigidity  of  the 
plate.   On  B1,  M„  and  V  are  related  to  M  .  M  and  M   by  [1], 

Mn  s  2  Mxyab  -  Mxa2  -  Myb2  ,    VR  s  Qn  -  A  ^  ,         (2, 


where 


Mns  ^Mx  "  V&b  +\y(a2-b2)  , 
Ti     9x  x   dy  xy      oy  y   dx  xy 


Here  and  in  the  remainder  of  the  paper  the  following  notation 
is  used:   subscripts  on  capital  letters  such  as  M,  Q  and  V  are 
Indices;  subscripts  on  lower  case  letters  such  as  u,  v,  w  and  <f> 
denote  partial  derivatives. 
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and  the  notation  df/ds  =  t  .  Vf  is  used  for  the  tangential 
derivative  of  f. 

Let  <J>(x,y)  be  the  unique  solution  of 

A<j)   =  -p    ,         for        x,y   in  R    ; 

(*) 

<J>   =   0,    on   C;         <J)      =   n.  V<f>   =  -V,    on  B*     . 

Here  A   is   the   two   dimensional   Laplacian.      The  boundary  value 
problem   (4)    is   called  Problem  <J>.      The  moment  potentials  u(x,y) 
and   v(x,y)    are   defined   by 


M      =   v     -  <j>    ,         M     =  u    -  h    ,        M        =    (v     +u    )/2    .  (5) 

x  y  y  xr'  xy         vxyy/ 


The  differential  equations  of  Southwell's  formulation  [2]  are 
then  given  by, 

Au  =  kV(-V.  u+2i>)  ,    k  =  (l-v)/(l+v)  ,         (6) 

where  u  is  the  vector  with  components  u  and  v  and  V  is  the  two 
dimensional  gradient  operator.   They  are  analogous  to  the 
displacement  equilibrium  equations  of  plane  stress  and  plane 
strain,  where  u  is  then  the  displacement  vector,  V<b  is  propor- 
tional  to  the  body  force  and  k  is  a  different  function  of  v. 
The  boundary  conditions  are, 

w  =  w  =  0  ,   on  C  ,  (7) 

Mn  =  M(s)  ,    Vn  =  V(s)  ,   on  B1  .  (8) 
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are  expressed    in   I 
(7)   implies  tha 
and  (5)   in  (1)   gives  the  re 


jbservlng   I 


v   v     =  u+  v=  <     . 


ux-v   vy  =  uy 


■  rsely,  if  u  satisfies  (9)  then  from  (5)  and  (1)  it  fol 
that  C  is  clamped  v;      a  rigid  rotation.   When  (5),  (2)  and 
(3)  are  inserted  in  (8),  then  (8)  is  reduced  to 

■  •  u  =  4>  -M  ,   |-  (n-  u  )  =  0  .  (10a, b) 

»v   «vS  OS   ~   ~s 


Here  the  notation  u  is  used  for  the  vector  with  components  u„ 

~s  s 

and  v  .   Equation  (10b)  is  integrated  with  respect  to  s  and  the 
s 

resulting  equation  and  (10a)  are  solved  for  u  .   Integrating  u 

gives 

s 

u  =  /   (<b-  M)t  ds  ,    on  B1  .  (11) 

~  J  ~ 


The  three  constants  of  Integration  are  set  equal  to  zero  since 
they  correspond  to  solutions  which  have  vanishing  bending  moments. 

Problem  P  is  defined  as  the  boundary  value  problem  which 
consists  of  determining  a  vector  u  that  satisfies  (6),  (9) 
and  (11).   The  solution  of  Problem  <{>  is  used  to  determine  the 
inhomogeneous  terms  in  Problem  P.   When  solutions  of  Problems  $ 
and  P  are  obtained  the  moments  are  determined  from  (5)«   The 
normal  displacement  can  then  be  evaluated  by  solving  (1)  for  the 
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second  derivatives  of  w  and  then  integrating  these  expressions 
using  (7).   Since  we  are  primarily  interested  in  the  moments, 
these  integrals  are  not  evaluated.   An  important  advantage  of 
Southwell's  formulation  is  that  the  complicated  higher  order 
derivative  conditions  on  B1  of  the  classical  formulation  are 
replaced  by  the  simple  Dirichlet  conditions  (11). 

At  the  corners  x  =  0,  and  y  =  yQ,  Jlf    the  boundary  condi- 
tions are  not  uniquely  prescribed;  see  (9)  and  (11).   Thus  we 
expect  a  singularity  in  the  solution  of  Problem  P  or  its  deriva- 
tives at  the  corners.   The  moments  may  be  discontinuous  and  even 
unbounded  at  these  points.   The  numerical  results  discussed  in 
Section  5  indicate  that  this  might  occur  for  the  specific  problems 
investigated. 


3.  The  Iterative  Method. 

An  approximate  solution  of  Problem  P  is  obtained  by  the 


following  iterative  method.   Starting  from  an  initial  estimate 

u  (x>y)  of  the  solution  for  fixed  data,  a  sequence  of  iterates 

1        2 
u  (x.y),  u  (x,y ),...,  is  defined  by  the  recursions 

Aun  =  kV(-V.  un-1+2cj>)    ,        x,y   in  R    ,  (la) 

^-vvj-1,        v^-u""1,        one,  (lb) 

s 
un  =  J       (<j>  -M)t   ds    ,         on  B*    ,  (lc) 

0 

un   =   e  un   +    (l-9)un_1    ,         in  R+B    .  (Id) 
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:-e  0   is  the  acceleration  parameter.   The  vector  u"  is  a 

ate.   The  function  <{>  in  (1)  is  a  solution  of 
Problem  <£ .   We  assume,  for  the  present,  that  it  is  known. 

Each  Iterate  in  (1)  is  determined  by  solving  the  boundary 
value  problem 

Aq  =  f(x,y)  ,    for  x,y  in  R  , 

q  =  h(y)  ,   on  C  ,    q  =  g(s)   for  x,y  on  B1  . 


(2) 


Here  f,  g  and  h  are  determined  by  <j>,  the  data  and  previous 
iterates.   An  approximate  solution  of  (2)  is  obtained  by  a 
difference  method  similar  to  the  one  employed  in  [7,8,10].   The 
region  R  is  covered  with  a  rectilinear  mesh  with  uniform  mesh 
spacing  6  in  the  x  and  y  directions.   In  this  paper  two  regions 
R  are  considered.   One  of  them  is  the  rectangle  and  the  other  is 
the  swept  back  plate  shown  in  Figure  1.   When  R  is  a  rectangle 
then  the  mesh  is  selected  so  that  mesh  lines  coincide  with  the 
edges  of  the  rectangle.   When  R  is  the  region  shown  in  Figure  1, 
then  x  =  0,  1  and  y  =  0,  1  are  taken  as  mesh  lines  and  each  mesh 
line  with  1  >  y  >   a  has  a  point  on  the  skew  edge.   For  more 
complicated  regions  it  may  not  always  be  possible  to  select  a 
uniform  mesh  such  that  mesh  points  coincide  with  the  boundary. 
Special  approximations  are  then  required  at  the  boundary. 

There  are  1+  1  and  J+  1  lines  in  the  x  and  y  directions 
so  that 

I  =  V6  '   J  =  v/6 '  (3) 
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When  R  is  a  rectangle,  I     and  H     are  the  lengths  of  the  sides. 

x      y 

When  R  is  the  region  shown  in  Figure  1  then  I     =   &     =1.   The 

x    y 

mesh  points  (x.,y.)  are  defined  by 


x.  =  i  6  ,    i  =  0, .  • .  ,1  , 


W 


y  •  =  j  6  ,   J  =  o,. . .  ,j  . 


The  set  of  mesh  points  which  are  in  R  is  called  Ro.   Likewise, 
the  set  of  mesh  points  on  B   is  called  Bo  and  the  set  of  mesh 
points  on  C  is  called  Co.   The  end  points  x  =  0,  y  =  y_  and 
x  =  0,  y  =  y,  are  assumed  to  be  in  Co.   A  function  defined  only 
at  the  points  of  the  mesh  is  called  a  mesh  function.   It  is 


assumed  that  at  each  point  x. ,  y.  of  Rc  +  Bc  + Cc  the  solution  of 

l  J  j     o    o    o 

(2),  q(x.,y.),  is  approximated  I 
solves  the  difference  equations 


(2),  q(x.,y.),  is  approximated  by  the  mesh  function  q. .  which 


ViJ   =  fiJ    '  (xi'yj)    e  R6    '  ^a) 


q±J  =  gu   ,         (x^yj)   e  B;  ,  (5b) 


D6qij   =  hij    '  (0  >7^    e  C6    '  ^5c) 


In  (5)  the  following  notation  is  used.   The  mesh  function  {g. .) 
is  defined  only  for  points  x. ,  y  .  in  Bc,   It  is  evaluated  by 
approximating  the  integrals  on  the  right  side  of  (1c)  by  Simpson's 
rule.   The  mesh  function  {h.  .)  is  defined  only  for  points  0,  y . 
in  Cg.   It  is  evaluated  by  approximating  the  derivatives  on  the 
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of  (lb)  by  a  centered     erence  formula.   When  R 
a  rectangle  the  mesh  f      n  (:'.  .]  La       d  by  appro- 
mating  the  de:  n  the       hand  side  of  (la)  by 
centered  difference  formulae.   If  R  is  the  region  in  Figure  1, 

■n  centered  difference  formulae  are  used  everywhere  except  in 
apprc  ed  second  derivative  along  the  diagonal 

e  Sc  which  consists  of  the  mesh  points  one  point  to  the  left 
of  the  skew  edge.   There  a  special  inward  difference  formula, 
with  the  same  truncation  error  as  the  centered  approximation,  is 
used  to  avoid  introducing  mesh  points  exterior  to  R.   Several 
different  inward  formulae  were  employed.   They  gave  only 
insignificant  changes  in  the  numerical  results.   The  difference 
operator  Dfi  in  (5c)  is  the  forward  three  point  difference 
approximation  of  0(6  )  of  d/5x.   In  (5a)  A&  is  the  nine  point 
difference  approximation  of  the  Laplacian  for  x.,y.  e  Re  -  Sc. 
On  Sfi  the  usual  five  point  Laplace  difference  operator  is  used. 
This  again  avoids  the  introduction  of  mesh  points  exterior  to  R 
that  are  required  if  the  nine  point  formula  is  employed  on  S&. 
The  solution  of  the  difference  equations  is  presumably  an 
accurate  approximation  of  the  solution  of  (2)  for  sufficiently 
small  6. 

The  difference  equations  (5)  are  a  system  of  linear 
algebraic  equations  for  the  mesh  function  (q.  .).   The  coefficient 
matrix  is  tridiagonal  with  respect  to  matrices.   The  system  is 
solved  by  a  direct  factorization  method  [9]«   Some  of  the  details 
of  the  application  of  the  factorization  method  are  discussed 
in  [7]. 


-10- 


Problem  <J),  (2A),    is  solved  numerically  in  the  same  way 
as  (2).   The  resulting  algebraic  equations  are  also  tridiagonal 
with  respect  to  matrices.   They  are  solved  by  the  factorization 
method  [9].   The  numerical  solution  of  Problem  <f>  is  used  to 
compute  the  inhomogeneous  terms  in  (1). 

When  the  iterations  (1)  have  converged,  the  moments  are 
computed  from  difference  approximations  of  (2.5). 


4.  Computing  Methods. 

A   sufficiently  small   value  of  the  mesh  spacing  5  was 
determined  by  performing  a  sequence  of  test  calculations  with 
successively  finer  meshes.   Mesh  sizes  of  5  =  l/lO,  l/l3,  l/26, 
l/40  and  6  =  l/52  were  used.   For  the  rectangle  the  mesh  size 
6  =  1/13  was  satisfactory.   The  test  calculations  for  the  swept 
back  plate  indicate  that  the  5  =  l/52  mesh  and  perhaps  a  finer 
mesh  is  required.   The  agreement  between  the  results  for  6  =  1/26 
and  6  =  l/52  is  good  for  mesh  points  away  from  the  corners 

x  =  y  =  0  and  x  =  0,  y  =  1.   Within  3  or  4  mesh  points  of  the 

** 
corners  the  results  are  significantly  different. 


The  computations  were  performed  on  the  CDC  6600  computer  at 
the  AEC  Computing  and  Applied  Mathematics  Center  of  the  Courant 
Institute  of  Mathematical  Sciences.   In  all  the  calculations 
v  =  0.32  was  used. 

Cf.  the  comments  in  Section  5« 
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The  following  numerical  convergence  criterion  was  employed: 

In  a  test  calculation  the  right  hand  side  of  (1)  was  decreased. 
This  Introduced  only  insignificant  changes  in  the  answers.  The 
condition  (1)  is  only  necessary  for  convergence. 

Starting  from  a  fixed  initial  iterate,  the  number  of  Itera- 
tions needed  to  satisfy  (1)  depends  on  the  value  of  the  accelera- 
tion parameter  9   in  (3.1).   The  optimal  value,  9        to  use  is  the 
one  that  minimizes  the  number  of  iterations  for  convergence. 
Estimates  for  9     were  obtained  from  a  series  of  test  calculations. 
For  the  problems  considered  9  =   0.6  was  a  satisfactory  value. 

The  number  of  iterations  for  convergence  varies  with  the 
data  of  the  problem  and  the  choice  of  9   and  6.   For  most  problems 
approximately  400  iterations  were  required.   For  a  mesh  with  676 
points  the  iterations  take  approximately  one  minute  of  CDC  6600 
time. 


5.  Applications  of  the  Method  and  Summary  of  the  Results. 

Some  of  the  numerical  results  will  be  compared  with  the 
exact  solution  of  the  infinite  strip  cantilever  0  <  x  <   I, 
|y|  <  00  ,  with  the  loads  independent  of  y.   Its  solution  which 
is  indicated  by  a  superscript  1  is  given  by, 
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I 

K   =f     t^-x)p(l)  d£-  (£-x)V  ,    MJ  =  V  Mj  ,    MJy  =  0  .     (1) 
x 

We  assume  that  M  =  0  in  (3.1c)  for  all  the  problems  studied. 

The  rectangular  plates  are  of  length  I   in  the  x  direction 
and  y  =  -y  =  -1.  If  p  and  V  are  symmetric  in  y  with  respect  to 
y  =  0,  then  M  and  M  are  symmetric  and  M   is  antisymmetric. 
The  square  cantilever,  I   =   2,  with  p  =  1  and  V  =  0  is  considered 
first.   Thus  B1  is  free.   The  results  shown  in  Figures  2  and  3 
are  obtained  with  6  =  l/l3.   The  maximum  moment  obtained  with 
this  mesh  is  M  (0,±2/l3)  £  -2.135.   The  graphs  in  Figures  2 
indicate  that  M1  is  an  accurate  approximation  for  x  >  0.3»   Near 
x  =  0  it  approximates  M  (x,0)  more  accurately  than  it  approximates 

-A. 

M  (x,-12/l3).   M   is  usually  much  smaller  that  the  other  moments 
x  xy 

and  M   =  0.   The  bending  moments  on  the  clamped  edge  C  are 
xy 

graphed  in  Figure  3.  The  approximation  (l)  is  accurate  in  the 
interior  of  C.  For  example,  it  underestimates  M  (0,0)  by  less 
than  6.5  percent. 

Graphs  of  M  on  the  clamped  edge  for  uniformly  loaded, 

.A. 

p  =  1,  V  =  0,  rectangular  plates  with  I   =  1,2( square)  K   are  shown 
in  Figure  4.   The  numerical  values  of  M  (0,0)  and  M  (0)  are 

J\.  A. 

summarized  in  Table  I. 
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Table  I.   Rectangular  Plates,  p  =  1,  V  =  0. 


6 

.0) 

-Mi(0) 

xx  ' 

°jb    error 

1 

1/20 

.514 

0.5 

-2.92 

2 

1/10 

2.135 

2.0 

-6.32 

2 

1/13 

2.135 

2.0 

-6.32 

2 

1/20 

2.136 

2.0 

-6.37 

k 

1/10 

8.968 

8.0 

-10.8 

For  the  wide  plate  £   =  1,   M  (0,y)  and  also  M  (0,y)  (not  shown) 

x  y 

have  three  stationary  points.   For  the  long  plate  £  =   4  and  the 
square  plate  £  =   2,  the  same  quantities  have  a  single  minimum. 
The  approximation  (1)  is  more  accurate  at  x  =  0  for  a  larger 
range  of  y  values  for  £  =   1   than  for  £  =   2  and  £  =   4.   In  fact, 
we  conjecture  that  as  £  —  0  the  solution  of  Problem  P  with  p  =  1 
and  V  =  o  approaches  (1)  in  the  interior  of  the  rectangle.   Near 
y  =  ±1  there  is  a  boundary  layer  in  which  the  solution  varies 
rapidly  from  (1)  so  as  to  satisfy  the  boundary  conditions  on 
y  =  ±1. 

The  values  of  M„(0,y)  for  £   =   2  and  6  =  1/10  and  6  =  1/20 
are  in  good  agreement  for  y  away  from  the  edges  y  =  ±1.   The 
values  of  MV(0,±1)  for  a  sequence  of  mesh  widths  are  listed  in 
Table  II. 
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Table  II.   Corner  Moments  for  the  Square 
Plate,  p  =  1,  V  =  0. 


6 

1/10 

1/13 

1/20 

Mx(0,±l) 

-.460 

-.269 

+  .037 

The  values  of  M  (0,±1)  Increase  as  6  decreases.   Because  of  the 
singularity  at  x  =  0,  y  =  tl,  it  is  not  clear  from  the  numerical 
results  whether  the  exact  value  of  M  (0,±1)  is  unbounded;  cf. 
the  discussion  in  [11]. 

Numerical  solutions  have  been  obtained  for  rectangular 
plates  that  are  deformed  by  a  uniform  shear  reaction  V  =  -1  on 
x  =  &f    |y|  *  1  and  tne  edges  y  =  ±1,  0  <_  x  <_  £   are  free.   The 
results  are  qualitatively  similar  to  those  of  the  uniformly  loaded 
cantilever  shown  in  Figures  2-4.   The  values  of  M  (0,0)  are  given 

-A. 

in  Table  III. 


Table  III.   Rectangular  Plates  with  p  =  0 
and  V  =  -lonx=,0. 


1 

6 

Mx(0,0) 

Mx(0) 

°/o    error 

1 

1/20 

1.057 

1 

-5.39 

2 

1/10 

2.220 

2 

-9.91 

2 

1/20 

2.219 

2 

-9.87 

4 

1/10 

4.644 

4 

-13.87 
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'."))  provides  an  underestimate.   The  accuracy  of  the 
estimate  increases  as  I   decreases.  The  corresponding  estimates 

.ire  more  accurate  for  the  uniformly  loaded 
rectangular  cantilevers;  see  Table  I. 

The  square  plate  with  B'  free  is  subjected  to  a  local 
surface  pressure  p  given  by 

H  ,    |x-x'|,  |y-y'|  <   6/2  , 

P  -  j  (2) 

^0  ,    |x-x'|,  |y-y'|  >   6/2  . 

Since  6  =  1/13,  the  total  force  produced  by  the  pressure  (2)  is 
F  =  6  =  I/169.   Representative  moments  obtained  from  the 
computations  with  x'  =  25/13,  y1  =  -11/13  are  shown  in  Figures  5. 
The  maximum  moment  is  M  (24/13,-1)  =  2.304  X  10"5.   Since  the 
solution  is  significantly  unsymmetric,  (1)  is  an  inaccurate 
approximation.   The  reduced  moment  m  is  defined  by, 


m  =  max  Mx(x,y)/  //  p  dxdy  .  (3) 

R        '    R 


It  relates  the  maximum  moment  to  the  total  applied  force.   Let 
m  and  m1  denote  the  values  of  m  respectively  for  the  uniform 
pressure  and  the  pressure  (2).   Then  the  stress  concentration 
factor  m'/ro  is  numerically  determined  as, 

m'/m°  =  .7295  .  (M 


The  swept  back  cantilever  plate  shown  in  Figure  1  is 

clamped  onx=0,  0  <  y  -  1  and  Bl  is  free.   Representative 

results  for  the  uniformly  loaded  plate  p  =  1  are  shown  in 

Figures  6.   The  values  of  M  (x,l/26)  obtained  with  the  meshes 

6  =  1/26  and  6  =  I/52  are  extremely  close  for  x  >  .2.   Near 

x  =  0  they  are  considerably  different.   This  is  presumably 

caused  by  the  singularity  at  x  =  y  =  0.   It  is  not  clear  whether 

M  (0,0)  is  bounded  as  &  -+  0.   The  values  of  M  (x,l/26)  obtained 
x  x 

with  the  6  =  1/26  mesh  are  in  qualitative  agreement  with  the 

corresponding  results  for  square  plate;  see  Figure  2a.   The 

values  of  M  (x,l/2)  obtained  with  the  two  meshes  are  within  the 

resolution  of  the  graph  in  Figure  6a.   Since  the  swept  back  plate 

is  geometrically  unsymmetric  M   (x,l/2)^  0.   However,  the 

xy 

twisting  moment  is  generally  much  smaller  than  the  other  moments 
and  is  not  shown  in  the  graphs. 

The  rapid  variation  of  M  (0,y)  near  y  =  0  as  6  decreases 

.A. 

is  presumably  a  result  of  the  singularity  at  the  corner.   The 
maximum  moment  is  |M  (0,6/52) |  £  .4o6.   The  maximum  values  of 
|M  (0,y)|  obtained  with  the  two  meshes  and  from  (1)  are  summarized 
in  Table  IV. 


Table  IV.   Maximum  Moment,  Uniformly 
Loaded  Swept  Plate. 


5 

max  |M  (0,y)| 
0<y<l 

MJ(0) 

9o  error 

1/26 

•  395 

.5 

27.2 

1/52 

.406 

•  5 

23.2 
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The  maximum  occurs  near  y  =  0  since  a  greater  force  is  exerted 

e  pressure  near  the  bottom  than  the  top  y  =  1  of  the  plate. 
The     :iite  strip  cantilever  therefore  inaccurately  predicts 
e  moments  at  x  =  0.   Let  m  denote  the  value  of  m   (3)  for 

5 

e  uniformly  loaded  swept  back  plate.   Then  the  numerically 
determined  stress  concentration  factor 


m°/m°  =  1.32  (5) 


occurs  when  the  swept  back  plate  is  required  to  support  the  same 

total  load  as  the  square  plate.   Since  the  maximum  moment  for 

uniformly  loaded  rectangular  plates  can  be  reasonably  approximated 

by  M  (0),  see  Table  I,  it  is  of  interest  to  determine  the 

variation  of  the  stress  concentration  with,  the  angle  of  sweep, 

the  aspect  ratio  and  the  quantity  a;  see  Figure  1.   We  have  not 

conducted  any  calculations  to  estimate  these  variations. 

1^(0, y)  for  the  swept  back  plate  with  the  pressure  (2)  with 

6  =  1/26  and  x'  =  23/26,  y'  =  1/13  is  shown  in  Figure  7.   The 

results  for  the  swept  back  plate  and  the  square  plate  with  the 

same  load  are  qualitatively  similar.   The  values  of  M  (0,y)  for 

the  square  plate  with  the  pressure  (2)  and  the  swept  back  plate 

with  a  uniform  pressure  that  gives  the  same  total  force  as  (2) 

are  also  shown  in  Figure  7-   This  gives  an  Indication  of  the 

stress  concentration  produced  by  localized  loads  and  sweep  back. 

Let  m'  denote  the  value  of  m  for  the  swept  back  plate  with  the 
s 

local  pressure.   Then 


-18- 


m's/m°Q   =   2.89  ,    m^/m'  =  1.313  .  (6) 

The  factors  in  (4),  (5)  and  (6)  suggest  that  the  combined  effects 
of  sweep  back  and  load  concentration  are  significantly  greater 
than  each  separate  effect. 

The  last  problem  discussed  is  the  swept  back  plate  with  the 
data 

V  =  0  ,    p(x,y)  =  y  .  (7) 

The  maximum  moment  on  the  clamped  edge  is  shown  in  Figure  8. 
The  ratio  M  /M  is  usually  smaller  for  the  load  (7)  than  for 

Xy/     X  J 

I 
the  uniformly  loaded  plate.   Let  m  denote  the  value  of  m  with 

s 

the  pressure  (7)«   Then 

m°/mf  =  1.75  .  (8) 

This  happens  because  the  larger  magnitudes  of  p  in  (7)  act  over 
a  smaller  area  and  hence  (7)  exerts  a  smaller  moment  at  x  =  0. 
Therefore  we  also  expect  M  (0,y)  to  be  more  symmetrical  about 
the  mid-line  y  =  1/2;  see  Figures  6b  and  8. 
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a=  1/13 


Figure  1 
Sketch  of  Swept  Back  Plate. 
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.5 


M(x,-l2/i3) 


.  | 1«-  y  =  - 1 1  / 1  3 


Figure  2a 

Variation  of  moments  on  the  line 
y  =  -12/13  for  the  square  plate  with 
a  uniform  pressure  p  =  1.   The  dashed 
line  is  the  moment  obtained  from  the 
infinite  strip  cantilever;  see  (5.1). 


Figure  2  b 

Variation  of  moments  on  the  center  line  y  =  0  for 
the  square  plate  with  a  uniform  pressure  p  =  1. 
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Figure  jj 

Variation  of  moments  on  the  clamped  edge,  x  =  0, 
the  square  plate  with  a  uniform  pressure  p  =  1. 
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Variation  of  the  moments  on  the  line  y  =  -11/13 
for  the  square  plate  with  a  localized  load  at 
x'  =  23/15,  y'  =  -11/13. 
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Figure  5b 

Variation  of  the  moments  on  the  center  line 
y  =  0  for  the  square  plate  with  the  local 
pressure. 
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Figure  5c 

Variation  of  the  moments  on  the  clamped 
edge  x  =  0  for  the  square  plate  with 
the  local  pressure. 
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Figure 


Variation  of  M  on  the  lines  y  =  l/26 

and  y  =  1/2  for  the  swept  back  plate 
with  a  uniform  pressure  p  =  1  and  the 
mesh  widths  6  =  1/26  and  6  =  1/52. 


Figure  6b 
Variation  of  M  (0,y)  on  the  clamped 

-A. 

edge  for  the  swept  back  plate  with 
p  =  1  and  the  mesh  widths  6  =  l/26 
and  6  =  1/52.   The  maximum  value  of 
M  (0,y)  for  6  =  1/26  is  -  .126  and 

occurs  at  y  =  2/13;  for  6  =  1/52  it 
is  -  .130  and  it  occurs  at  y  =  J>/26. 
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Figure  7 

Curve  1  is  M  (0,y)  for  the  swept  back 

plate  with  a  local  pressure  (5.2) 
where  6  =  1/26  and  x!  =  23/26, 

Curve  2  shows  the  corre- 
sponding M  (0,y)  for  the  square  plate 

with  the  same  pressure.   Curve  3  is 
M„(0,y)  for  the  swept  back  plate  with 

-A. 

a  uniform  pressure  which  gives  the  same 
total  force  as  the  local  pressure  . 
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Figure  8 
The  variation  of  M  (0,y)  on  the  clamped  edge  for 

-A. 

the  swept  back  plate  with  the  linear  pressure 
distribution  (5-7).   The  graph  of  M  (0,y)  is 

similar  to  the  one  for  M  (0,y)  but  with  small 

magnitude.   The  maximum  value  of  M  (0,y )=- .02477 
and  occurs  at  y  =  2/13.  y 
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This  report  was  prepared  as  an  account  of 
Oovernment  sponsored  work.   Neither  the 
United  States,  nor  the  Commission,  nor  any 
person  acting  on  behalf  of  the  Commission: 
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the  information  contained  in  this  report, 
or  that  the  use  of  any  Information, 
apparatus,  method,  or  process  disclosed 
In  this  report  may  not  infringe  privately 
owned  rights;  or 

B.  Assumes  any  liabilities  with  respect  to 
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